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We present systematic calculations of the properties of 11~ isomers in neutron-deficient 184 ~ 198 pb 
and iss-^Po. These states 

are based on the 7r(/ig/ 2 «i3/2)K' r =ii- configuration. They are calculated 
in the framework of the Hartree-Fock-Bogoliubov method with a Skyrme interaction and density- 
dependent pairing force. The energies and deformations of the 11~ states are compared to those of 
the intruder + states in both Pb and Po isotopes. In the most neutron-deficient Po isotopes, the 
calculations predict, below a weakly oblate II - state, another oblate 11~ state which is even more 
deformed than the intruder oblate 11~ state in their Pb isotones. The energies and quadrupole 
moments of the 11~ isomers, corresponding to a weakly oblate nuclear shape, are in fair agreement 
with the available experimental data. 

(3JT)' PACS numbers: 21.60.Jz,21.10.Dr,21.10.Re,21.10.Ky,27.80.+w 
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The neutron-deficient nuclei around the Z = 82 shell closure constitute one of the most famous examples of shape 
coexistence 0,0 ■ I n Pb isotopes, oblate, prolate and superdeformed configurations coexisting with a spherical ground 
state were predicted theoretically (see [3j and references therein) and identified experimentally 0, 3 more than ten 
years ago. By now, low-lying + states, which are a characteristic feature of shape coexistence, are known from 202 Pb 
down to 184 Pb H, IE S 13 • The interpretation of these low-lying + states relies on proton excitations across the 
Z = 82 spherical shell gap. In a shell model picture 0, EJ, the energy of these multi-particle- multi-hole (mp-mh) 
configurations is lowered by the quadrupole-quadrupole and pairing interactions. In mean-field models |lllll2l |. these 
Cf") ' structures are described by a state with an oblate (prolate) quadrupole deformation, whose low energy is due to the 
occupation of deformed intruder orbits, such as the down-sloping 7r/ig/ 2 -orbital for an oblate state. 

High-spin isomers constructed on these + states are of particular interest since the measurement of their static 
electromagnetic moments provides direct information on both their single-particle structure and their deformation. 
G ' Several isomeric states are known in Pb isotopes. In particular, 11~ levels have been observed from 188 Pb to 198 Pb 
(see, e.g. Ref. ^3(). They are not known in heavier Pb, most probably because they are not yrast. These states are 
interpreted as a two-quasiparticle (2qp) excitation based on the ft (hg^Ha/ 2) ic* =u- single-particle configuration as 
confirmed in 196 Pb by the measurement of its o-factor Experimental information on deformation is extracted from 
the collective bands based on the excited + |M Fil Il4l Il5| and II - states Recent data on static quadrupole 
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moments for 11~ states in 194 > 196 p D [T8HI{^ | support the mean-field interpretation of these states as weakly deformed 
• »-h oblate configurations. 

The experimental excitation energies of the 11~ states and of the 0^" states in Pb isotopes are shown in the left 
part of Fig. 1. Both exhibit a similar behavior. The steady decrease of energy with decreasing neutron number 
indicates that, in a shell model picture, the residual quadrupole and pairing interactions, or, in mean-field models, 
the deformation effects, are maximal at the middle of the neutron sub-shell. 

Similar coexisting structures based on intruder configurations are known experimentally above and below the Z = 82 
shell closure, in the neutron-deficient Hg, Pt and Po isotopes 0, 0] . The energy systematics of the low-energy excited 
+ states in neutron-deficient 20 °- 196 p o 0,[2(| strongly resembles that in Pb. The observation of bands based on 
these states seems to indicate a strong mixing between different shapes. The 11~ isomers have been identified from 
the semi-magic 210 Po down to 194 Po |l3|. In contrast to the Pb case, their excitation energy is nearly constant as a 
function of the neutron number. 

Previous theoretical studies of these isotopes have been mainly devoted to the description of the + intruder states, 
in particular in Pb. A large variety of models has been used: a schematic spherical shell model [a. |9j. l21l a nd its 
symmetry-based boson truncations [22|. the Strutinsk y appro ach with a Woods-Saxon potential [ill Il2l l23j. and 
several variants of self-consistent mean- field methods |24l EEl I2H l27j . A comprehensive study of the low-energy 
spectroscopy and ground-state properties of light Po isotopes has been presented within several schematic approaches 
in Ref. |2g (see review 0] for more references). On the other hand, the only systematic study of isomers in these 
neutron-deficient nuclei is based on cranked Nilsson-Strutinsky type calculations with a schematic Woods-Saxon 
potential and is devoted to 186 - 200 Pb [ill PbSf. The quadrupole moments of ll - intruder states in Pb and Po down 
to 194 Pb and 198 Po have also been analyzed recently within the particle-core coupling model [p|. The deformation 
properties of the ll - state in 196 Pb have been discussed in the framework of the tilted-axis cranking model [lg- The 
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magnitude of the ll - nuclear deformation is indeed a basic ingredient in this model for the description of magnetic 
rotational bands in both even- A and odd-A neutron-deficient Pb isotopes 0, |2|j. However, up to now, no fully 
microscopic investigation of these isomeric configurations has been performed. 

In this letter, we present the results of the first fully self-consistent description of the properties of 11~ intruder 
isomers in 184 - 198 p D an d in 188 -200p o k a g ec j on the Hartree-Fock-Bogoliubov (HFB) method. Our aim is to analyze 
the nature of these states and to study the evolution of their deformation as a function of neutron number. The 
simultaneous consideration of Pb and Po isotopes will allow also to determine the effect of the Z = 82 spherical shell 
closure on the deformation trends. Since several Po isotopes have extremely shallow energy surfaces as a function 
of the axial quadrupole moment, we have also performed selected configuration-mixing calculations to validate the 
construction of isomers on an oblate configuration. 

All our calculations have been performed with the same set of nucleon-nucleon effective interactions. In the mean- 
field channel, we use the Skyrme parameterization Sly4 |30| and in the pairing channel, a density-dependent zero-range 
force mainly active on the nuclear surface as introduced in Ref . [3lJ . The HFB equations have been solved on a three- 
dimensional cubic mesh for reflection symmetric nuclear shapes with the method presented in Refs. [32l and 
include an approximate particle- number projection with the Lipkin-Nogami prescription |34j . The isomeric states 
are constructed fully self-consistently as 2qp excitations as in Ref. 35], selecting the quasiparticles whose dominant 
components in the HF basis are the 9/2~ and 13/2 + particle states. 

In Pb isotopes, we have taken for the energy and the deformation of the first excited + states those of the oblate 
minimum of the deformation energy curves. Configuration- mixing calculations pil l25l l26l |27| have shown that this 
approximation is well justified, even when the potential well of the oblate minimum is very shallow. The topology of 
the energy curves for Po isotopes is more complex, with several shallow minima on both prolate and oblate sides. The 
properties of the first excited + state have therefore been determined by performing configuration- mixing calculations 
with respect to the axial quadrupole collective coordinate. To avoid the ambiguity resulting from the fact that the 
states obtained from the configuration-mixing calculations do not have the same mean-particle number, we have 
included in the configuration mixing a projection on the correct proton and neutron numbers. 

This set of methods has already been applied with success to several other problems, such as the description of 
neutron-rich nuclei (3^, H|| , of superdeformed rotational bands [3^, 13]} and of heavy nuclei (3^ . 

Deformation energy curves 

The deformation energy curves obtained for 194 - 184 pb and 198 - 188 p isotopes are shown in Figs. 2 and 3, respectively. 
These curves have been obtained by HF+BCS calculations with a constraint on the axial mass quadrupole moment 
q = <720 = (2z 2 — x 2 — y 2 ). Our formalism for configuration mixing calculations is indeed up to now limited to the 
mixing of HF+BCS wave functions and does not permit to use HFB ones. For this reason, we shall base our discussion 
of the deformation properties on the HF+BCS energy curves. We have checked that very similar curves are obtained 
within the constrained HFB method and that the trend in the 0^" excitation energy as a function of the neutron 
number is the same. In order to make a connection with simpler models, we have labelled orbitals by approximate 
Nilsson indices. 

For all Pb isotopes (Fig. 2), the ground state corresponds to the spherical minimum. With decreasing neutron 
number, two structures appear, first as shoulders, then as secondary minima on the oblate (starting from A = 192) 
and prolate sides (starting from A = 188). The oblate minimum occurs around q = — 10 b which is the deformation 
corresponding to the crossing between the down-sloping 9/2[514] mean-field orbital and the up-going orbital l/2[440]. 
On the prolate side, the minimum at q ~ 20 b corresponds to the deformed shell gap due to the 1/2 [530], 1/2 [541] and 
3/2[532] intruder orbitals and the rising of the l/2[400], 3/2[402] and ll/2[505] orbitals. For 186 Pb, the excitation 
energy of the prolate minimum significantly overestimates the experimental data. This discrepancy is irrelevant for 
the description of the isomeric states which are constructed on the oblate minimum. The energy of this prolate state 
is very sensitive to the details of the mean-field method. Several other mean-field calculations, based either on the 
HF+BCS method and other Skyrme interactions 0i i or on the HFB method with the Gogny force [2!| , predict 
a prolate minimum lower in energy than the oblate one. 

As we have explained in the introduction, for the Pb isotopes, we have approximated the energy of the oblate + 
state by the energy of the oblate minima. These energies are compared with the experimental data in the left part of 
Fig. 1. The decrease of the HF+BCS energies with decreasing neutron number is in good agreement with the data, 
although the experimental energies are overestimated. This overestimation is even larger by approximately 300 keV 
in HFB calculations. 

The deformation energy curves for Po isotopes (Fig. 3) exhibit a complicated coexistence of several shallow minima. 
To check whether these minima are stable against triaxial deformations, we have performed selected calculations 
including constraints on the triaxial quadrupole moment. These calculations show that the symmetric minima 
corresponding to q ~ ±8 b (j3 ~ ±0.1) obtained for all Po isotopes are separated by a very low barrier of about 100 
keV in the 7 plane. They are therefore probably unstable against 7 vibrations. The shoulder that appears in 198 Po 
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for a deformation around —15 b becomes a deep well for lighter isotopes which is the ground state for 192 > 190 p o . In 
188 Po, the oblate minimum is nearly degenerate in energy with a prolate minimum at q ~ 19 b. The Strutinsky 
calculations of Refs. ^jH^,^ agree qualitatively with our results. It is clear from this discussion that a study of the 
low lying levels of Po isotopes must take into account the vibrational motion with respect to the axial quadrupole mode. 

Properties of low-lying oblate + states: configuration-mixing calculations for Po isotopes 
We have performed configuration- mixing calculations with respect to the axial quadrupole moment for the Po isotopes 
within the generator coordinate method (GCM), as was presented in Refs. [3jj|40j. The basis for the calculations is 
a set of wave functions Q(q) obtained by HF+BCS calculations with a constraint on the axial quadrupole moment. 
These wave functions $(<?) are projected on the correct number of protons and neutrons. Collective wave functions are 
then constructed as a linear superposition of the projected non-orthogonal basis functions. The weights corresponding 
to the different quadrupole moments are determined by diagonalizing a discretized version of the Hill and Wheeler 
equation |4lll42|. This procedure provides as many eigenstates as discretization points of the quadrupole moment. The 
lowest eigenvalue corresponds to the ground state with a wave function spread over the axial quadrupole moment. 
As we shall see, depending on the structure of this wave function, it represents either a vibration within a well- 
defined well or a mixing of mean-field states with different shapes. In the same way, the other eigenvalues of the 
Hill and Wheeler equation represent excited states. Depending upon their shape, they can be classified into three 
main categories. Some states are mainly located in a well that is clearly separated from the ground state (as it is the 
case for a super-deformed well). Other states can be [3 vibrations with respect to a lower state in the same potential 
well. Finally, some states can represent a mixing of states corresponding to different shapes. We have not included 
triaxial deformations in these configuration-mixing calculations. Some Po isotopes being soft against the 7 degree 
of freedom, these deformations may affect some of our results. However, their inclusion would considerably increase 
the computation time, and since our main purpose is to calculate the 11 _ isomers, we have limited the configuration 
mixing to axial shapes. 

The energies of the first few states obtained by the diagonalization of the Hill and Wheeler equation are shown for 
selected Po isotopes in the upper part of Fig. 4, together with the deformation energy curves. The triangles representing 
the energy of the configuration- mixed states are positioned at their average quadrupole moment. The amplitudes of 
the collective eigenf unctions (which are related to the mixing coefficients by an integral transformation (4^) are plotted 
in the lower part of the figure for the three lowest states obtained for each isotope. The wave functions of these states 
are spread over a large range of quadrupole moments. In 20 °. 198 : 196 p O; the ground state and the first excited state 
are a mixing of an oblate and a prolate configuration, resulting in an average small oblate deformation. This mean 
quadrupole moment is close to the one of the first shallow oblate minimum of the deformation energy curve. In the 
three isotopes, the second excited + state has an average quadrupole moment q ~ — 10 b, corresponding to a weakly 
oblate shape. While it has only a slightly dominant weight in the ground-state wave function of 200 < 198 . 196 p O] this 
oblate configuration becomes the ground state in 194 > 192 > 190 p Oj with still a spreading over a large range of quadrupole 
moments. In 190 Po, it is nearly degenerate with a predominantly prolate configuration. These results agree well 
with the interpretation of the experimental data that the ground states of 20 °- 196 p o isotopes are mixed but nearly 
spherical, the ground state of 194 Po is mixed, with a large oblate component and that of 192 Po is predominantly 
oblate p ,l44 | . Experiments also suggest a strong mixing between oblate and prolate shapes in the 190 Po ground state 

(Ellil) 

In the right part of Fig. 1, we compare the three first eigenstates of the configuration- mixing calculation with the 
available experimental data. The states that have the largest oblate mean deformation are indicated by filled triangles. 
The experimental excited + states are known in 200 > 198 < 196 p o [2(j and are interpreted as 4p-2h excitations in the 
spherical shell model (see Refs. and references therein) or as states with an oblate deformation in mean-field 

models. The experimental trend is well reproduced by the most oblate states, although it is not the first excited state. 
Note however that the experimental energies are intermediate between those of the two first excited states. This may 
indicate that the interaction that we use leads to a too large mixing between the deformed configurations. Such a 
result could be modified by the inclusion of triaxial deformations, but in a way which can hardly be anticipated. 

From the deformation energy curves for the Po isotopes, the 11~ states can be constructed on either of the two 
oblate minima. The configuration-mixing analysis clearly favors the slightly oblate minimum for the heavier isotopes. 
For the lighter isotopes ( 194_188 Po), the mean deformation of the configuration-mixed ground state is intermediate 
between the two. For these nuclei, we shall therefore explore whether an II - isomer can be constructed for both 
deformation regions. 

Properties of the ll - isomers 
The 11~ states in the Pb and Po isotopes have been constructed as a 2qp excitation on an HFB vacuum. As in Ref. [35| . 
the 2qp state is calculated fully self-consistently, which means that effects such as core polarisation due to the time- 
odd terms of the interaction and the conservation of the number of particles or changes of deformation are taken into 
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account. Since the single-particle basis, which is formed by the eigenstates of the Hartree-Fock Hamiltonian, does not 
diagonalize the qp basis, the qp excitations cannot be defined exactly by single-particle states. We have always chosen 
the 2qp wave functions which have the largest overlaps with the lowest 9/2~ and 13/2+ single-particle states. In 
practice, these overlaps are always larger than 0.95. As in previous calculations of lqp and 2qp excited states H3,|35|, 
the isomeric state has been obtained axial in all cases, although triaxial deformations are allowed by the symmetries 
of our computer code. 

The experimental and theoretical excitation energies of the 11~ states are compared in Fig. 1 for the Pb and 
Po isotopes. In the Pb isotopes, the energy of the 11~ states decreases with decreasing neutron number, both in 
experiment and theory. The experimental trends are quite well reproduced by our calculation, although the excitation 
energy is overestimated by 600 to 800 keV. This overestimation is most probably related to a better description in 
the HFB formalism of the vacuum (corresponding to the ground state) than of qp excitations. Our calculation shows 
that the excitation energy of the 11~ state is the lowest in 190 > 188 pb and increases in 186 < 184 pb for which there are no 
experimental data yet. 

For the Po isotopes, the ll - states are known down to 194 Po. The isomer has been identified in 192 Po 44], but 
its energy could not be obtained from the data. In 200 - 196 p O) a minimum corresponding to an 11~ state could be 
constructed only for small oblate deformations. The theoretical excitation energies are in good agreement with the 
data. In the case of 194 - 188 p 0; 2qp states have been obtained for both oblate deformations. The isomers corresponding 
to smaller (larger) absolute values of the quadrupole moment are shown by filled (empty) diamonds in Fig. 1 (right 
part). The excitation energies of both solutions are close to each other. These states have to be further mixed by a 
configuration-mixing calculation. One should not expect a significant change in the excitation energy of the lowest 
mixed eigenstate, the ground state being lowered in a similar way by configuration mixing. 

The energies of the 11~ states are nearly constant as a function of the neutron number and the isomers occur at a 
lower energy than in Pb. The 11~ energies do not follow the steep decreasing trend as the intruder + energies. This 
suggests that the structure of the 11~ isomers is rather similar for all Z = 84 isotopes, at least down to 196 Po. In 
lighter Po isotopes, the 11~ states built on the most oblate minimum become lower in energy than the weakly oblate 
ones, and configuration- mixing calculations would be necessary to get more accurate predictions on the 11~ energy 
and deformation. 

In the terminology of phenomenological models, the constant energy of the II - states in heavier Po isotopes, 
indicates that the self-consistency effects due to the 2qp excitation do not affect strongly the core. This is not the 
case for the Pb isotopes, in which the Z = 82 magic number is broken by the 2qp excitation, inducing an additional 
core deformation, as also demonstrated by the trend in the intrinsic quadrupole moments. 

The calculated intrinsic charge quadrupole moments of the 11~ states in Pb and Po are plotted in Fig. 5. For 
comparison, there are also given two available data points for 194 . 196 pb extracted from the measured spectroscopic 
quadrupole moments in the assumption of the axial symmetry and K = II 18] . The theoretical results are in 
good agreement with the experimental values. The deformation parameters (/32, (3a) defined as in Ref. |46| vary for 
Pb isotopes from (-0.15,-0.014) for 198 Pb to (-0.19,-0.018) for 184 Pb. The quadrupole moments of the intruder 
isomers in the Pb- isotopes are significantly larger than those of the 11~ states corresponding to the weakly oblate 
shape in the heaviest Po isotopes. This reflects the influence of core polarisation due to the breaking of the Z — 82 
core in the Pb isotopes. 

In lighter Po isotopes, two candidates for an 11~ isomer have to be mixed by a configuration-mixing calculation. 
It is difficult to predict the quadrupole moment of the mixed state which could have an intermediate value between 
those of the two mean-field states. From these results, one can infer a smooth change from a slightly deformed oblate 
shape of ll - states in the heaviest Po to a more deformed oblate shape in the lightest. Experimental observation of 
the rotational bands above these isomers or a direct measurement of the spectroscopic quadrupole moments would 
be helpful to clarify the situation. For the unambiguous ll - states in the heavier Po isotopes, the deformation 
parameters (j3 2 , (3 4 ) vary from (-0.114, -0.007) for 200 Po to (-0.117, -0.007) for 196 Po. 

In summary, we have, for the first time, applied the HFB method with a Skyrme interaction and a realistic pairing 
force to the study of the 11~ isomers in a scries of neutron-deficient Pb and Po nuclei. The isomers have been 
constructed self-consistently as 2qp excitations 7r(/i9/2-*i3/2+)n- • Their main properties are in good agreement with 
the available experimental data. We obtain that the energy of 11~ states in Pb isotopes reaches minimum around 
N = 106 — 108 and increases in lighter isotopes for which there are no data yet. For Po isotopes, the ll - energy is 
predicted to remain roughly constant as a function of neutron number. The calculations support the assumption of 
oblate deformation of these nuclei in their II - states, with a steady increase towards neutron-deficient isotopes. The 
deformation is more pronounced in Pb than in Po, but in lighter Po isotopes the deformation of the ll - states might 
become of the same order or even larger than that in the Pb isotopes. We have also studied the properties of the 
low-lying + states in Po isotopes within the configuration-mixing approach on the basis of HF+BCS wave functions. 
These calculations show that, in all cases, the lowest + states do not have a well defined deformation, and oblate, 
spherical and prolate shapes are strongly mixed, both in the ground and the excited + states. 
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FIG. 1: Theoretical and experimental energies of the 11" and states in 184 " 198 p D (left) and in 188 - 200 p o (right) (see Ref. [3 
for the references). The experimental 0^ in 184 Pb is not indicated, since its nature is not well established Q. For 184 Pb, the 
oblate + state from configuration-mixing calculations is shown for comparison. For the theoretical + states in Po isotopes, 
the first three eigenvalues obtained in the configuration-mixing calculations are represented by triangles. The filled triangles 
correspond to the states with the largest oblate mean deformation. The empty diamonds denote the 11~ states based on the 
second oblate minima for the cases when convergence was achieved. 





FIG. 2: Deformation energy curves for Pb obtained from the constrained HF+BCS calculations. 
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FIG. 3: Deformation energy curves for 198 ~ 188 p obtained from the constrained HF+BCS calculations. 




FIG. 4: Energies of the GCM eigenstates, together with the HF+BCS deformation energy curves for 198 - 190 p o (upper part). 
The amplitude of the corresponding collective GCM eigenfunctions for three lowest states enumerated by k are shown in the 
lower part of the figure. 
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FIG. 5: Comparison between the charge quadrupole moment of the 11~ states in Pb (triangles) and Po (diamonds). The empty 
diamonds correspond to the 11~ states constructed on the second oblate minima in light Po isotopes. The intrinsic quadrupole 
moment as extracted from the experiment [Tsl Il9| for 196 > 194 Pb are shown by circles with error bars. 



